Charm production asymmetries from heavy-quark 
recombination 



^ \ Thomas Mehen f 



Department of Physics, Duke University, Durham NC 27708, USA 
PACS numbers: 13.85.Ni, 14.40.Lb, 14.20.Lq, 12.39.Hg 



O 
O 

(N 
£ 

Abstract. Charm asymmetries in fixed-target hadroproduction experiments are 
>> \ ', sensitive to power corrections to the QCD factorization theorem for heavy quark 

production. A power correction called heavy-quark recombination has recently been 
proposed to explain these asymmetries. In heavy-quark recombination, a light quark 
^ ' or antiquark participates in a hard scattering which produces a charm-anticharm quark 

pair. The light quark or antiquark emerges from the scattering with small momentum 
in the rest frame of the charm quark, and together they hadronize into a charm particle. 
' The cross section for this process can be calculated within perturbative QCD up to 

. an overall normalization. Heavy-quark recombination explains the observed D meson 

| and A c asymmetries with a minimal set of universal nonperturbative parameters. 

P-i! In this talk I will discuss asymmetries in the production of charm particles in fixed- 
Ph! target hadroproduction experiments ^ El IH Ej • These asymmetries are much larger 

than predicted by perturbative QCD calculations, so they are a sensitive probe of 
nonperturbative aspects of heavy particle production. Asymmetries in the production 
of light particles such as pions and kaons have also been observed, but the production 
of particles containing heavy quarks is under better theoretical control. When the mass 
of the quark is heavier than Aqcd, as is the case for charm, perturbative QCD can be 
applied even when the particle is produced without large transverse momentum. In 
addition, nonperturbative effects can be organized in an expansion in AqcdA?V 

The differential cross-section for the production of a charm particle, H, in the 
collision of two hadrons, A and B, is believed to factorize in the following manner jH] 

da[A + B -» H + X] = (1) 
fi/A ® fj/B ® da[ij -> cc + X] ® D C ^ H + ... . 

hi 

Here i,j denote partons, fi/A and fj/B are parton distribution functions, D c ^h is 
the fragmentation function for c hadronizing into H and da[ij — > cc + X] is a 
perturbatively calculable short-distance cross section. The ellipsis represents corrections 
to the factorized form of the cross section that are suppressed by Aqcd/^c, or possibly 
Aqcd/p_l if P± > m c . 
f mehen@phy.duke.edu 
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The leading contribution to the factorization theorem in Eq. (JJJ) predicts that charm 
particles and their antiparticles will be produced symmetrically. The leading order 
partonic processes, gg — > cc and qq — > cc, produce charm quarks and anticharm quarks 
with identical kinematic distributions. The fragmentation functions D c _,h and D-_jj 
are identical due to the charge conjugation invariance of the strong interactions. No 
asymmetry between H and H is generated at leading order in perturbation theory. 
Next-to-leading order perturbative corrections |H1 El E2] can generate asymmetries 
but these are quite small. The asymmetry, defined by 

^ = «fj, (2) 
a[H\ + a[H \ 

is never larger than a few percent. 

The asymmetries observed in experiments can be much larger El El El El • One 
well-known phenomenon is the "leading particle effect" in fixed-target hadroproduction. 
Charm hadrons sharing a valence parton with the beam hadron are produced in greater 
numbers in the forward direction of the beam than particles that do not share a valence 
quark with the beam. For instance, experiments with a n~ beam incident on a nuclear 
target observe many more D~ and D° than D + and 75° in the forward direction of the 
7r~ beam. a[-D~] rises from nearly zero in the central region, xf ~ 0, to 0.7 at the 
highest xp measured. In the forward region of the n~ beam almost six D~ are produced 
for every D + . Large asymmetries are also observed in charm baryon production. In pN 
collisions a[A+] m 1 for all ij? > 0.2. 

Smaller asymmetries are also observed even when there is no leading particle 
effect. Examples include fixed-target photoproduction jTTJ H21 ES], where the beam 
has no valence quark quantum numbers, the production of D s mesons in fixed-target 
hadroproduction experiments with tt or p beams, and A c production in experiments 
with 7r~ beams. In the last case, there is no asymmetry from the leading particle effect 
because the A+ shares a d quark with the n~ while the A~ shares a u with the tt~. These 
asymmetries are smaller than those observed when there is a leading particle effect but 
still larger than perturbative QCD predictions. 

Traditionally the asymmetries have been explained by nonperturbative models of 
hadronization. A commonly used model is the Lung string fragmentation model ^3] 
which can be implemented using PYTHIA [T^]. (Another model of string fragmentation 
can be found in Ref. [H3)- The asymmetry is generated by the "beam drag effect" ^7] in 
which the charm quark binds to the remnants of the incident hadron via the formation 
of a color string. These models can be tuned to fit the data but this can require 
an unusually large charm quark mass and large intrinsic transverse momentum for the 
partons in the incoming hadrons pQ . The PYTHIA Monte Carlo with default parameters 
rarely predicts the asymmetries correctly £Q and in the case of A c asymmetries in 
nN collisions |3] gets the sign of the asymmetry wrong. Another approach is the 
recombination model first introduced in Ref. [TBI (f° r recen t analyses using similar 
methods, see Refs. 01201121]). In this model the charm quarks coalesce with spectator 
partons in the beam hadrons whose momentum distribution is determined by double 



Charm production asymmetries from heavy-quark recombination 



Figure 1. Example of a diagram for cq recombination into a D meson. Single 
lines are light quarks, double lines are heavy quarks and the shaded blob represents 
hadronization into a state including the D meson. 

parton distributions. Finally there are models that rely on the existence of intrinsic 
charm in the incident hadron wavefunction to generate the asymmetry [221 EH EI]- These 
models are sensitive to a number of poorly determined nonperturbative functions such 
as distributions of partons in the remnant and functions that parametrize recombination 
probabilities. 

Since the asymmetries observed in nature are much larger than what is predicted 
by perturbative calculations, they are a direct probe of the power corrections to the 
QCD factorization theorem. Recent work has demonstrated that an O(AqcdA^c) 
power correction called heavy-quark recombination provides a simple explanation of 
charm meson and baryon asymmetries in photo- and hadroproduction experiments 
[23 EE1 EH EE| This approach differs from previous nonperturbative models in that 
the asymmetry is generated in the short-distance process so cross sections are calculable 
up to an overall normalization set by a few universal nonperturbative parameters. For 
D mesons, the dominant contribution comes from a process in which a light antiquark 
participates in a hard scattering process that produces a charm- ant icharm quark pair. 
The light antiquark emerges from the hard scattering with momentum of O(Aqcd) in 
the rest frame of the charm quark, then the light antiquark and charm quark hadronize 
into a final state that includes a D meson. The process is depicted in Fig. [0 This cq 
recombination process gives a contribution to the D meson cross section of the form 

da[qg — > D + X] = 2J da[qg — > cq(n) + c] p[cq{n) — > D] . (3) 

n 

In this formula, da[qg — > cq(n) + c] is a short- distance cross section for producing a cq 
with quantum numbers denoted by n and p[cq(n) — > D] parametrizes the hadronization 
of the cq into a state that includes a D meson. 

The NLO correction to the fragmentation cross section in Eq. (0) includes a 
subprocess qg — > qcc which is similar to the short-distance part of the mechanism 
depicted in Fig. ^ However, in the NLO calculation the c and q hadronize 
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Figure 2. Example of a diagram for cq recombination, the leading mechanism for 
charm baryon production via heavy-quark recombination. 

independently, so this correction does not account for the possiblility that the cq can 
bind nonperturbatively and hadronize into the D. For this reason, the heavy-quark 
recombination mechanism is an O(AqcdA^c) correction to Eq. (0) rather than part of 
the NLO correction to the fragmentation contribution. 

Because the light antiquark is massless it is natural to expect the heavy-quark 
recombination contribution to be a convolution of a short- distance cross section with a 
distribution function that depends on the fraction of the light-cone momentum carried 
by the light quark. However, to lowest order in Aqcd/^c only the leading moment 
of such a distribution contributes. Because the cross section is inclusive, the final 
state may include other soft quanta besides the D meson. Soft gluons emitted in 
the hadronization process can change the total angular momentum and color quantum 
numbers of the cq produced in the short-distance process. Therefore the color and 
angular momentum quantum numbers of the cq can be different from the D meson. 
Amplitudes for production of cq in L > partial waves are suppressed by additional 
powers of Aqcd/^c relative to S- waves and can be neglected. It is also possible for 
the light antiquark flavor quantum number of the D to be different than that of the 
cq due to light quark-antiquark pair production in the hadronization process. However, 
light quark-antiquark pair production is suppressed in the large- N c limit of QCD so this 
effect is a subleading correction to D meson production. This leaves four parameters 
for heavy-quark recombination into D + mesons: 

p 1 = p[cd( l s! ) 1) )^D + }, p 1 = p[cd(*s[ 1) )^D+], (4) 
p 8 = p[cdCS®) D + ], p 8 = P [cd( 3 S[ 8) ) -> D+}. 

These parameters scale with the heavy quark mass as Aq CD / m c and explicit expressions 
in terms of nonperturbative QCD matrix elements can be found in Ref. |29|. Analogous 
parameters for D° and D~ mesons are obtained by using isospin symmetry and charge 
conjugation invariance, while parameters for D* + states are related to those in Eq. (JH) 
by heavy-quark spin symmetry. 
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To form a baryon from cq recombination requires the production of at least 
two light quark-antiquark pairs which is suppressed by 1/N%, so cq recombination 
should contribute mostly to D meson production. The process which gives the 
leading contribution to charm baryon production is cq recombination which is shown 
in Fig. 121 This is the same as cq recombination except a light quark participates in 
the recombination process instead of a light antiquark. A light quark-antiquark pair 
must be produced for the cq to hadronize into a state that has a charm baryon or 
charm meson, so there is a 1/N C suppression in either case. Thus, cq recombination 
is the dominant heavy-quark recombination contribution to baryon production but is 
subleading for meson production. In this talk I will compare heavy-quark recombination 
predictions with data on A c asymmetries. The heavy-quark recombination contribution 
to the A+ production cross section is 



The 7] parameters are similar to the p parameters described earlier. There are two 
possible color states, 3 and 6, and two possible spin states contributing at this order, 
for a total of four parameters for A+ production: 



Isospin symmetry requires r)[cu(n) — > A+] = rj[cd(n) — > A+]. These parameters also 
scale with the heavy quark mass as AqcdA?V 

Techniques for calculating the short-distance cross sections da[qg — > cq{n) + c] and 
da[qg — > cq(n) + c] as well as explicit expressions can be found in Refs. [23 EH]- All 
these cross sections are strongly peaked in the forward direction of the initial light quark 
or antiquark. Define 9 to be the angle between the light quark (or antiquark) in the 
initial state and the cq (or cq) in the center-of-momentum frame. The ratio of the 
recombination cross section to the dominant leading order fragmentation cross section, 
gg — *> cc, at 9 = it/ 2 is 



where s is the parton center-of-mass energy squared. When the cq or cq emerges at 
right angles to the incoming partons, heavy-quark recombination is suppressed relative 
to fragmentation by a kinematic factor of m 2 c /s. It is also suppressed when the charm 
quark is produced in the backward direction, 9 = it: 



da[qg -> A+] = d °[l9 -»■ cg(n) + c] r][cq(n) -> A+] . 



(5) 



773 = vM'Sn - A+], 773 = r/[cu( 3 Sj 3) ) -> A+], 
V6 = vicuCS^) - A+], fk = v[cuCs[ 6) ) - A+]. 



(6) 





da[qg^>cq( 3 sP) + c] 
da[gg -> cc] 



da[qg -> cq^S?) + c] 
da[gg -> cc] 




(8) 



da\qg -> cq^S^) + c] 
da[gg -> cc] 



da[qg -> cq^S^) + c] 
da[gg -> cc] 




(9) 
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Note that the backward suppression is greater when the cq or cq is produced in a 1 Sq 
state. In the forward direction, there is no kinematic suppression: 



da[qg — > cq{n) + c] 



da[qg —> cq{n) + c] 
0=0 da l99 -> cc] 



(10) 

6»=0 



da[gg -> cc] 

Heavy-quark recombination preferentially produces charm mesons and baryons in the 
forward direction of the incident light quark (or antiquark). The cross section is 
larger for charm particles that share a valence quark with one of the colliding hadrons 
because the structure functions of the valence quarks is largest. Thus, the heavy-quark 
recombination mechanism provides a natural explanation of the leading particle effect. 

Another mechanism by which charm hadrons can be produced is by ordinary 
fragmentation of charm quarks produced in cq or cq recombination. This process, 
called "opposite side recombination" , gives the following contributions to charm hadron 
production: 

da[qg -> H + X] = £ da[qg -»■ cq{n) + c] p[cq(n) ^D}® D C ^ H , (11) 

n,D 

da[qg -> H + X] = ^ d ^9 ~* cq(n) + c] rj[cq(n) -»■ B ] <g> D C ^ H . (12) 

Here B is a charm antibaryon. The opposite side recombination mechanism can generate 
asymmetries even when there is no leading particle effect. 

Fig. El compares calculations of the D ± asymmetry using the heavy-quark 
recombination mechanism with data from the E791 experiment in which a 500 GeV 
7r~ beam is incident on a nuclear target |2j. In this analysis the feeddown from D* 
mesons is included but feeddown from excited D meson states is neglected. The charm 
quark mass is set equal to 1.5 GeV and the renormalization and factorization scales are 



•p\ + m 2 c . The parton distributions are GRV 98 LO [20] for the nucleon and GRV-P 
LO [SI] for the pion. The strong coupling constant used is the one-loop expression for 
a s with 4 active flavors and Aqcd = 200 MeV. Only the LO fragmentation diagrams 
are included in the calculation. The effect of NLO perturbative corrections can be taken 
into account by multiplying the LO cross section by a K factor of about 2 (THj. The 
fragmentation functions are (^-functions times fragmentation probabilities taken from 
Ref. At fixed-target energies the single particle inclusive distributions da/dxF and 
da/dp\ are better reproduced with 5-function fragmentation functions rather than, for 
example, Petersen fragmentation functions jlOj . 

In this calculation of the D meson asymmetry the opposite side cq recombination 
contribution of Eq. (|T2~j) is not included, so the result is independent of the r] parameters. 
The p parameters were determined by a global analysis of all measurements of D meson 
asymmetries in fixed-target hadroproduction experiments [22] • Fits with one, two and 
three parameters were performed. (A four-parameter fit did not yield significantly better 
results than the three-parameter fit.) The results of all three fits are shown in Fig. |21 
Note that even a fit with one parameter, p 1 = 0.14 and all other p parameters set to 
zero, does a good job of describing the data in the forward region where heavy-quark 
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Figure 3. a[D~] vs. xp and p\ for a 500 GcV tt~ beam on a nuclear target One-, 
two- and three-parameter fits are solid, dashed and dot-dashed curves, respectively. 

recombination is most important. Inclusion of other parameters is necessary to obtain 
agreement in the central region xp ~ 0. 

The universality of the p parameters is tested by looking at asymmetries of other 
D mesons as well as asymmetries in experiments with different beams. Fig. H] compares 
the predictions of the heavy-quark recombination mechanism for a[.D* _ ] with data from 
the E791 experiment j2]. Fig. shows data and a calculation of a[-D~] for the WA89 
experiment which has a 340 GeV X~ beam incident on a nuclear target These 
calculations use the same parameter sets as in Fig. Eland the agreement with experiment 
is excellent. 

Finally, the cq recombination mechanism for charmed baryons is tested by 
comparing to measurements of «[A+] in experiments with 500 GeV 7r~ beams from 
E791 |3] and 540 GeV p beams from SELEX The results are shown in Figs. Eland 
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Figure 4. a[D* ] vs. xp for a 500 GcV ir beam 0. Solid, dashed and dot-dashed 
curves are the same as Fig. [3] 
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Figure 5. a[D s ] vs. xp for a 340 GeV £ beam Solid, dashed and dot-dashed 
curves are the same as Fig. [21 

[71 In the analysis of baryon asymmetries, the single nonvanishing p parameter is chosen 
to be consistent with the one-parameter fit to D meson asymmetries. Setting all the 
rj parameters to zero gives the asymmetry shown by the dotted lines in Figs. El and UJ 
which is generated entirely by opposite side recombination. Though this is adequate 
for the 7T~ beam data, the p beam data clearly requires an additional mechanism for 
generating baryon asymmetries. A one parameter fit with 773 = 0.22 is shown by the 
solid lines in Figs. El and UJ The results of the calculation are in good agreement with 
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Figure 6. a[A+] vs. xf for a 500 GeV 7r~ beam 0]. The solid curve is the best 
single-parameter fit with rj 3 = 0.22, while the dotted curve is in the absence of cq 
recombination. 
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Figure 7. a[A+] vs. xf for a 540 GeV p beam Fit parameters for the solid 
and dotted curves are the same as Fig. The horizontal line at a = 1 is the physical 
upper bound. 

the baryon asymmetry measured in both experiments. 

In this talk, I have shown that the heavy-quark recombination mechanism provides 
a natural mechanism for generating the observed charm meson and baryon asymmetries. 
There are four cq recombination parameters and four cq recombination parameters 
which are universal. The heavy-quark recombination mechanism correctly accounts 
for observed D meson and A c asymmetries. It is especially encouraging that the same 
set of p parameters can be used to correctly describe D, D* and D s asymmetries in two 
different experiments with different beams. A c asymmetries in both tt~ and p beams 
are well described by a fit with a single p parameter that is consistent with D meson 
data and a single 77 parameter. In the future, a more detailed analysis of charm hadron 
asymmetries which thoroughly tests the universality of the p and rj parameters should 
be performed. Once these universal parameters are determined, predictions for heavy- 
quark recombination cross sections for charm and bottom hadron production in a large 
variety of nuclear and particle physics experiments will be possible. 
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